Measurement of Giant Dipole Resonance width at low temperature: A 

new experimental perspective 

S. Mukhopadhyay , Deepak Pandit Surajit PaP, Srijit Bhattacharya^, A. De'^, 
S. Bhattacharya^, C. Bhattacharya'^, K. Banerjee'^, S. Kundu^, T. K. Rana^, 
G. Mukherjee^, R. Pandey^, M. GohiP, H. Pai^, J. K. Meena^ and S. R. Banerjee'^'* 

'^Variable Energy Cyclotron Centre, 1/AF-Bidhannagar, Kolkata-700064, India 
^Department of Physics, Barasat Govt. College, Barasat, N 2^ Pgs, Kolkata - 700124, India 
'^Department of Physics, Raniganj Girls' College, Raniganj - 713358, India 



Abstract 

The systematic evolution of tiie giant dipole resonance (GDR) width in the temperature region of 0.9 ~ 1.4 MeV has 
been measured experimentally for ^^^Sb using alpha induced fusion reaction and employing the LAMBDA high energy 
photon spectrometer. The temperatures have been precisely determined by simultaneously extracting the vital level 
density parameter from the neutron evaporation spectrum and the angular momentum from gamma multiplicity filter 
using a realistic approach. The systematic trend of the data seems to disagree with the thermal shape fluctuation 
model (TSFM) . The model predicts the gradual increase of GDR width from its ground state value for T > MeV 
whereas the measured GDR widths appear to remain constant at the ground state value till T ^ 1 MeV and increase 
thereafter indicating towards a failure of the adiabatic assumption of the model at low temperature. 

Key words: Low temperature GDR width; Adiabatic thermal shape fluctuation model; BaF2 detectors. 
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In recent years, great interests are being shown 
to understand the correct description of damping 
mechanisms contributing to the Giant Dipole Res- 
onance (GDR) width at low temperatures (T) [1- 
3]. The damping of such giant collective vibration 
inside the nuclear medium occurs either due to es- 
cape of resonance energy by means of particle or 
photon emission (escape width) or due to its redis- 
tribution in other degrees of freedom within the sys- 
tem (spreading width) [4]. In medium and heavy 
nuclei, it turns out that the escape width only ac- 
count for a small fraction and the major contribu- 
tion of the large resonance width comes from the 
spreading width [5,6]. The general trend of the res- 
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onance width, as deduced from the lorentzian fit to 
the cross-section data, has been found to be smallest 
for the closed shell nuclei and larger for the nuclei be- 
tween shells [4]. However, it needs to be mentioned 
that for deformed nuclei, the width was obtained by 
fitting one single Lorentzian. Such fits never resulted 
in a systematic mass dependence of the width. Re- 
cently, an empirical formula has been derived for 
the spreading width by separating the deformation 
induced widening from the spreading effect and re- 
quiring the integrated Lorentzian curves to fulfil the 
dipole sum rule [7]. The relation has been found to 
hold good for the widths of the different GDR com- 
ponents corresponding to the three axes of a de- 
formed nucleus in general [7,8]. The inclusion of the 
deformation in describing the apparent GDR width 
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is also supported by the recent experimental [9] and 
theoretical [10] development in the description of the 
nuclear ground states. 

The GDR, built on excited state, is an important 
experimental tool since it couples directly to the 
nuclear shape and the investigation of its strength 
distribution gives a direct access to the nuclear 
deformation. Owing to this property, it has been 
applied to study Jacobi shape transition [11,12] 
and hyper-deformation [12] in alpha cluster nuclei. 
This triaxiality is caused by J-drivcn deformation 
in high J and T regime and could be measured ex- 
perimentally since the shapes are characterized by 
very large deformation (/3 > 0.6). The experimental 
results have been substantiated with a theoretical 
calculation based on the thermal shape fluctua- 
tion model [TSFM] [13], which takes into account 
the J-driven deformation and T-driven shape fluc- 
tuation. For a triaxial non-rotating nucleus, the 
GDR strength function is a superposition of three 
Lorentzians that correspond to the vibration of 
the nucleus along each of the semi-axes [14]. The 
resonance energy corresponding to the each axis is 
obtained using the Hill- Wheeler parametrization 

Ek = Eq exp (7 — 27rA;/3) while 

the widths are calculated applying the power law 
Tfc = ro{Ek/Eoy-'^ [15]. Eo and To are the param- 
eters for a spherical nucleus with mass A. Further, 
at very high angular frequencies, these three GDR 
components split (the ones perpendicular to the 
spin axis) due to Coriolis eff'ect as the GDR vibra- 
tions in a nucleus couple with its rotation when 
viewed from a non-rotating frame giving rise to five 
GDR components altogether [16,11,12]. Finally, 
the GDR cross-section is calculated by taking into 
account the large amplitude thermal fluctuations 
using a Boltzmann probability e~^^^'^'^^'^ with the 
volume clement (3* ii'm{3j)d(3d'-f, where F is the 
free energy [15]. However, for small deformations, 
experimentally it is not possible to measure the 
shape of the nucleus since thermal shape fluctuation 
smears out the associated splitting of the strength 
function resulting in an overall broadening of the 
distribution. Thiis, only the; apparcint GDR widths 
are measured from the experiment using a statisti- 
cal model analysis and compared with the TSFM, 
which also provides the apparent width of the GDR 
(in turn the shape of the nucleus). A systematic 
study of the thermal fluctuation model revealed the 
existence of a universal scaling law for the apparent 
width of the GDR for all T, J and A [17,18]. 



The apparent width of the GDR, built on the 
excited states, has been found to increase mono- 
tonically (^T^/^) [19] beyond T > 1.5 MeV. One 
should expect a gradual increase in the apparent 
GDR width from its ground state value (T = MeV) 
with the increase in temperature as predicted by 
TSFM. However, the temperature region below 1.5 
MeV has rarely been investigated to verify if such 
a behavior is really true. In Sn and nearby nuclei ( 
A ~ 120), mostly investigated so far, only a single 
apparent GDR width measurement exists for T < 
1.2 MeV which lies well below the TSFM prediction 
[1] . On the other hand, the phonon damping model 
(PDM) [20] which considers the couphng of the GDR 
phonon to particle-particle and hole-hole configu- 
rations as the mechanism for the increase of GDR 
width, without any need of T-driven shape fluctua- 
tions, attributes this suppression to thermal pairing 
which contributc;s oven beyond 1 MeV. These two 
models clearly disagree with one another at temper- 
atures below 1.5 MeV highlighting the importance of 
microscopic efifects responsible for this unusual phe- 
nomenon. In order to address these issues and to test 
the validity of the theoretical models, a systematic 
comparison between experiment and theory over a 
range of temperature for several nuclei is required. 

Experimentally, the measurement of GDR width 
at low temperature is very challenging due to the 
difficulties in achieving low excitation energy. Tra- 
ditional heavy ion fusion reactions are limited to 
higher temperature due to the presence of Coulomb 
barrier in the entrance channel and are always as- 
sociated with broad J distributions. Inelastic scat- 
tering [1.21,22] has been used as an alternative ap- 
proach with the advantage that the angular mo- 
mentum transfer will be relatively low, but, the ex- 
citation energy windows are uncertain to about at 
least 10 MeV and hence, the estimated temperatures 
are less precise. Due to these reasons, very few and 
widely separated (~ 0.25 MeV) data points with 
large error bars are available. In the present work, 
alpha induced fusion reactions with precise experi- 
mental techniques has been used to investigate the 
low temperature region. In these reactions, the de- 
scription of excitation energies from where the GDR 
photons decay will be more precise and the asso- 
ciated maximum angular momentum for ^ ~ 119 
mass region will be rather small (~ 207i). 

In this Letter, we present the first systematic and 
precise experimental study of angular momentum 
gated apparent GDR width in the unexplored low 
temperature region (0.9 - 1.4 MeV) for ^^^Sb us- 
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ing fusion reaction with alpha particles. The nuclear 
level density (NLD) parameter is a critical input for 
statistical model calculation and is also important 
for proper estimation of the nuclear temperature [4] . 
The NLD parameter has been measured experimen- 
tally from the neutron evaporation spectrum while 
the bremsstrahlung contribution has been estimated 
using the forward/backward 7-ray anisotropy [23]. 

The experiment was performed at the Variable 
Energy Cyclotron Centre (VECC), Kolkata using 
accelerated alpha beams from the K130 Cyclotron. 
A self-supporting 1 mg/cm^ thick target of ^^^In 
(99% purity) was bombarded with beams of '*He. 
Three different beam energies of 30, 35 and 42 MeV 
were used to form the compound nucleus ^^^Sb at 
the excitation energies of 31.4, 36.2 and 43.0 MeV, 
respectively. The LAMBDA high energy photon 
spectrometer [24] (98 large BaF2 detectors arranged 
in two blocks of 7x7 each) was used to measure the 
high energy gamma rays (> 4 MeV) at the angles 
of 55, 90 and 125 degrees with respect to the beam 
axis. The detector arrays were positioned at a dis- 
tance of 50 cm from the target. Since, the GDR pa- 
rameters depend on both the excitation energy and 
the angular momentum populated, it is important 
to separate the two effects in order to understand 
their individual contribution. Hence, along with the 
LAMBDA spectrometer, a 50-element low energy 7 
multiplicity filter [25] was used (in coincidence with 
the high energy 7-rays) to estimate the angular mo- 
mentum populated in the compound nucleus in an 
event- by- event mode as well as to get a fast start 
trigger for time-of-flight (TOF) measurements. The 
filter was split into two blocks of 25 detectors each 
which were placed on top and bottom of a specially 
designed scattering chamber at a distance of 5 cm 
from the target in a staggered castle type geometry. 
The TOF technique was used to discriminate the 
neutrons from the high energy 7-rays. The pulse 
shape discrimination (PSD) technique was adopted 
to reject the pile-up events in the individual detec- 
tor elements by measuring the charge deposition 
over two integrating time intervals (30 ns and 2 fis) 
[24]. The neutron evaporation spectra were mea- 
sured using seven liquid scintillator (BC501A, 5" 
diameter and 7" long) based neutron time of flight 
detectors [26] in coincidence with the multiplicity 
filter. The neutron detectors were placed at the an- 
gles of 30, 45, 75, 90, 105, 120 and 150 degrees with 
respect to the beam direction and at a distance of 
150 cm from the target. The time resolution of the 
neutron detectors was typically about 1.2 ns which 
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Fig. 1. [Color online] (top panel) The experimental gamma 
ray energy spectrum (blue circles) at 35 MeV is compared 
with the CASCADE prediction with bremsstrahlung com- 
ponent (black continuous line). The individual CASCADE 
(green dashed) and bremsstrahlung (red dotted) components 
are also shown, (bottom panel) The experimental ai coeffi- 
cient (blue circles) is compared with the exponential fit (red 
continuous line) with the same Eq parameter as used for the 
bremsstrahlung shape. 



gives an energy resolution of about 15% at 1 MeV 
for the present setup. 

In alpha induced fusion reactions at these ener- 
gies, the populated angular momenta are quite low 
to have any effect on the widening of the GDR width. 
The angular momentum dependent increase of the 
apparent GDR width starts showing up above a spin 
of 35^ for A - 120 mass region ( J - 0.6^^/^) [17]. 
However, it is still important to measure the angu- 
lar momentum to properly estimate the tempera- 
ture of the compound nucleus. The measured fold 
distribution from the multiplicity filter was mapped 
onto the angular momentum space using a Monte 
Carlo GEANT3 [27] simulation. The procedure is 
described in detail in ref [25] . The extracted angular 
momentum values for three different incident ener- 
gies are listed in Table 1. 

The energy of the evaporated neutrons has been 
measured using TOF technique whereas the neutron 
gamma discrimination was achieved by pulse shape 
discrimination (PSD) and TOF. The neutron TOF 
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spectra were converted to neutron energy spectra 
using the prompt gamma peaks in the TOF spectra 
as the time reference. The efficiency correction for 
the neutron detectors were done using the Monte 
Carlo Computer code NEFF [26]. The evaporated 
neutron energy spectra, after transformation from 
laboratory frame to centre of mass frame, were com- 
pared with CASCADE [28] calculation using Chi- 
square minimization technique in the energy range 
of 2 - 8 MeV for the determination of the level density 
parameter (a = A /a). The extracted inverse level 
density parameters (a) are listed in Table 1 with an 
uncertainty of ± 0.4 MeV. 

The high energy 7-ray spectra were generated 
in the offline analysis using the cluster summing 
technique [24] in which each detector element was 
required to satisfy the prompt time gate and pulse 
shape discrimination gate. The measured high en- 
ergy 7-ray spectra at 90° were compared with a 
modified version of the statistical model code CAS- 
CADE [28] along with a bremsstrahlung compo- 
nent. The non statistical contributions to the ex- 
perimental 7 spectra arising due to bremsstrahlung 
emission were parametrized using the relcition fx hrem 
= k/[C+exp(E^/Eo)] [23]. The centre of mass 7- 
ray angular distributions were assumed to have the 
form cr(6')=Ao[l + aiPi(cos(6')) -I- a2P2(cos(6'))] 
as the emission of gamma rays is dominated by 
electric dipole radiation [23]. The ai coefficient 
should be zero for statistical emission, however, 
it is non zero for higher gamma energies due to 
bremsstrahlung emission. The slope parameter (Eq) 
of the bremsstrahlung shape was extracted by si- 
multaneously fitting the ai coefficient using the ex- 
ponential function with the same slope parameter. 
The bremsstrahlung component as well as the ai 
coefficient for 35 MeV incident energy is shown in 
Fig 1. The extracted value of the slope parameters 
are consistent with the systematics Eq — l.l[(Eia6 
- Vc)/Ap]*'-^^, where Eiab, Vc and Ap are the beam 
energy, coulomb barrier and the projectile mass re- 
spectively [29] . The Coulomb barrier in the studied 
reaction is 15.1 MeV. In the CASCADE calcula- 
tion, the level density prescription of Ignatyuk [30] 
has been taken with the asymptotic level density 
parameter as extracted from the corresponding 
neutron evaporation spectrum. The simulated spin 
distributions deduced from the experimental multi- 
plicity distributions were used as inputs for different 
folds. The predictions from the CASCADE calcu- 
lations and the bremsstrahlung contributions were 
convoluted with the detector response and com- 
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Fig. 2. [Color online] Divided plots of the measured 7-spectra 
(blue circles) and the best fitted CASCADE calculations (red 
continuous lines) for different folds (F) at incident energies 
of 30, 35 and 42 MeV. 

pared with the experimental gamma-ray spectra 
for different folds. The best fit was obtained using 
a minimization technique in the region of 8 - 20 
MeV. In order to highlight the GDR region, both 
the data and calculated spectra were linearized by 
dividing with a statistical spectrum assuming con- 
stant El strength (Fig 2). The average temperature 
was estimated using the relation (T} = [(E* - E^ot 
- ^gdr) /o]^^'^ where E* is the excitation energy. 
Erot is the energy bound in the rotation at the aver- 
age J corresponding to a particular fold. The GDR 
centroid energies {Ejgdr) were found to be constant 
at around 15 MeV. The extracted parameters for 
different incident energies are listed in Table 1. At 
this point, it needs to be mentioned that the nu- 
clear deformation was not included in our analysis 
and we report on the extraction of the apparent 
GDR widths and compare them with the TSFM, 
which also provides the apparent width of the GDR 
including all shape fluctuations. 

The apparent GDR widths measured in the low 
temperature range of 0.9 - 1.4 MeV in the present 
study are shown in Fig 3 along with other mea- 
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Tabic 1 

GDR width and bremsstrahlung parameters at various folds for 30, 35 and 42 MeV beam energies. 



MeV 



30 



35 



42 



Fold 



2 

> 3 



2 

3 
> 4 



2 

3 
> 4 



h 



13.5 ± 4.7 
18.0 ± 5.2 



12.2 ± 4.7 
15.1 ± 4.8 
20.1 ± 5.2 



16.2 ± 4.8 
19.1 ± 4.9 
23.5 ± 5.3 



Erot 

MeV 



2.5 
4.3 



2.1 
3.1 
5.3 



3.5 
4.8 
7.1 



(T) 
MeV 



0.98 
0.93 



1.17 
1.13 
1.03 



1.36 
1.30 
1 .23 



Apparent GDR 
width (MeV) 



4.7 ± 0.3 
4.6 ± 0.3 



5.1 ± 0.3 

4.8 ± 0.3 

3.9 ± 0.3 



5.8 ± 0.3 
5.3 ± 0.3 
5.0 ± 0.3 



Eo 
MeV 



2.8 
3.0 
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Fig. 3. [Color online] Plot of apparent GDR width with tem- 
perature. The red filled circles are the values for ^^®Sb de- 
duced from the present work. For comparison, data for ^^"Sn 
from previous works are shown. The black downward trian- 
gle is from ref [1]. The green upward triangles and blue open 
circles are from ref [17] while the brown open squares are 
from ref [23] . The pink continuous line represents the TSFM 
calculation [17] while the black dotted line corresponds to the 
phonon damping model calculation [20]. The green dashed 
line is the apparent ground state width of ^^"Sb (discussed 
in the text). 

surements done earlier for ^^'^Sn [1,17,23]. The con- 
tinuous line represents the adiabatic thermal shape 
fluctuation calculation [17] at low spin. It is evident 
that the temperature dependence of the apparent 
GDR width determined from this experiment dif- 
fers substantially from the adiabatic thermal shape 
calculation at low temperature. In ^ ~ 120 mass 
region, where shell effects are small, the apparent 
GDR width is expected to increase with tempera- 
ture from its ground state value in a manner con- 
sistent with the properties of hot liquid drop [17]. 
In contrast, the systematic experimental data show 
the apparent GDR width to be constant till T ^ 1 
MeV and increases thereafter. The extracted appar- 



ent GDR widths at temperatures T < 1 MeV, match 
pretty well with the apparent ground state width 
of "^Sb (4.5 MeV, dashed line in Fig 3) as calcu- 
lated using the spreading width parametrization [7] 
r = O.OSEqI)^ for the small ground state deforma- 
tion (/3 — -0.12) [31]. The discrepancy between the 
experimental data and TSFM indicate the failure 
of the model in the present form in describing the 
evolution of the apparent GDR width with temper- 
ature below 1.5 MeV. A similar suppression of the 
apparent width compared to TSFM was observed, 
in the mass region ~ 117, at still lower temperature 
(0.68 MeV) by measuring the high gamma rays from 
the hot fission fragments produced in ^^^Cf cold fis- 
sion [3]. At these low temperatures, several micro- 
scopic effects might play significant roles and should 
be incorporated properly for a better explanation of 
the experimental data. But, even after incorporat- 
ing these corrections, the situation does not improve 
[17]. It is therefore needed to have a re-look into the 
formulation of the model with the incorporation of 
any microscopic effects that may be responsible for 
such deviation. 

We remark here that the present work establishes 
the fact that the apparent GDR width remains 
constant at ground state value till T ^ 1 MeV and 
increases subsequently thereafter as the apparent 
widths are measured below and above T=l MeV. 
This was not done in ref[l] and the conclusion was 
based on assumption, as only a single apparent 
GDR width (with large error) at T=l MeV was 
reported. Moreover, the data are measured simul- 
taneously with neutron evaporation spectra to put 
a constraint on the vital level density parameter 
in statistical calculation. The bremsstrahlung com- 
ponent in the high energy region is also extracted 
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experimentally instead of keeping it as free parame- 
ter. As a result, the present study provides a precise 
experimental systematics of the apparent GDR 
width in low T regime which provides an important 
test to the theoretical model predictions. 

The microscopic phonon damping model (PDM) 
[20] (the dotted curve in Fig 3), though not used 
widely, better explains the trend of the data at this 
low temperature region. The model calculates the 
GDR width and the strength function directly in 
the laboratory frame without any need for an ex- 
plicit inclusion of thermal fluctiiation of shapes. It 
has been shown that the thermal pairing effect plays 
an important role in lowering the GDR width at 
T <2 McV [20]. However, the PDM does not have a 
built-in angular momentum dependence of the GDR 
width at finite temperatures and may, therefore, be 
used only to describe the temperature dependence 
at zero spin. Nevertheless, it is interesting to find 
that the PDM, which does not take into account 
the deformation of the nucleus, better explains the 
data compared to TSFM that emphasizes on the 
inclusion of amplitude fluctuations in shape (defor- 
mation). This is in complete contrast to the GDR 
built on the ground states as the apparent widths 
are successfully described by including the nuclear 
ground state deformation. Both the models describe 
the data and show similar behavior for ^^"Sn (Fig 3) 
at higher temperatures ( T > 1.5 MeV) but at lower 
temperatures show drastically diff'erent behavior in- 
dicating, perhaps, the presence of more dominant 
quantal effects. This exciting result opens up a new 
question regarding the basic understanding of the 
increase of the apparent GDR width in the complex 
quantal nuclear many body system at low tempera- 
ture and is left for further theoretical insight. 

In summary, we have presented a systematic ex- 
perimental study of the apparent GDR width in the 
temperature region of 0.9 ~ 1.4 MeV in ^^^Sb using 
fusion reaction with alpha particle to explore the un- 
charted region of the evolution of GDR width with 
temperature. The apparent GDR widths deduced 
from the experiment are inconsistent with the pre- 
dictions of thermal shape fluctuation model [TSFM] . 
In fact, the GDR width appears to be constant at 
its ground state value till T ~ 1 MeV and increases 
thereafter whereas the TSFM predicts a gradual in- 
crease of GDR width from its ground state value for 
T > MeV. The discrepancy between the exper- 
imental data and TSFM can be conjectured as an 
artifact of quantal effects of damping mechanism at 
low temperature. 
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